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Abstract. The properties of matrix valued polynomials generated by the scalar- 
type Rodrigues' formulas are analyzed. A general representation of these polyno- 
mials is found in terms of products of simple differential operators. The recur- 
rence relations, leading coefficients, completeness are established, as well as, in 
the commutative case, the second order equations for which these polynomials are 
eigenfunctions and the corresponding eigenvalues, and ladder operators. 

The conjecture of Duran and Griinbaum that if the weights are self-adjoint and 
positive scmidcfinite then they are necessarily of scalar type is proved for Q = x and 
Q = x 2 — 1 in dimension two, and for any dimension under genericity assumptions. 

Commutative classes of quasi-orthogonal polynomials are found, which satisfy 
all the properties usually associated to orthogonal polynomials. 



1. Introduction 

Orthogonal matrix valued polynomials have been the subject of substantial in- 
vestigation in recent years, due not only to their intrinsic interest, but also to their 
applications. A beautiful introduction in the topic and a wealth of references can be 
found in [I], [2], [3], and in the recent overview [I]. 

The present paper investigates the properties of matrix valued polynomials gener- 
ated by the scalar-type Rodrigues' formulas. These have also been considered in [5] 
under the assumption that the weight is self-adjoint and positive semidefinite. The 
present paper investigates the properties of resulting polynomials under no assump- 
tions on the "weight"; in fact such polynomials appear in classification problems of 
differential equations [6]. 

It is shown that the scalar-type Rodrigues' formulas can be reformulated as a 
product of simple differential operators, Proposition!]], and their commutation rela- 
tions yield a very simple derivation of the recurrence relations, leading coefficients, 
completeness, and, in the commutative case, the second order equations for which 
these polynomials are eigenfunctions, and the corresponding eigenvalues, §2.41 



2 



RODICA D. COSTIN 



The problem of finding in which cases the "weights" are indeed weights, in the sense 
that they are self-adjoint and positive semidefinite, is investigated in the present pa- 
per for Q(x) = x and Q{x) = x 2 — 1. It is found that under appropriate nonresonance 
conditions (satisfied by generic matrices) the conjecture of Duran and Griinbaum in 
[S] is true: self-adjoint weights reduce to scalar ones (Proposition^ in the sense 
that for some S we have W(x) = SA(x)S^ 1 with A(x) diagonal. In two dimen- 
sions the same is true for most of the resonant cases (Proposition^. There exists 
one (resonant) case when there are self-adjoint weights not reducible to scalar ones 
(Proposition!!]); interestingly, all the generated polynomials belong to a commuta- 
tive algebra of matrices. However, these weights are positive semidefinite only when 
reducible to scalar ones; this completes the proof of Duran- Griinbaum conjecture in 
dimension two. 

The generated polynomials in this last case form a quasi-orthogonal family which 
satisfies all the usual properties associated to orthogonal polynomials, and they can 
be generalized to higher dimensions, §2.51 

2. Main Results 

2.1. The product representation. In this paper V denotes a finitely-dimensional 
complex vector space, M. denotes the matrices C(V,V), A4[x] are the .M-valued 
polynomials, and let Ai(x) denote the set of .M-valued functions. 
Let Q(x) be a polynomial degree at most two: 

(1) Q(x) = ax 2 + rx + 5 

(real or complex-valued) and let Li,L 2 G M. with L\ satisfying the (nonresonance) 
assumption: 

(2) L\ + ka is invertible in M. for all k — 1, 2, . . . 

Let W(x) be an M.- valued function satisfying the Pearson equation 

(3) Q(x)W(x)- 1 W'(x) = xLx + L 2 

Definition Let P n (x) be the .M-valued function defined by the Rodrigues formula 

(4) P n (x) = Wix)' 1 — [Q(x)»W(x)] 

Note that P = I and P\(x) = {2a + Li)x + r + L 2 . Proposition[TJ below shows 
that in fact P n {x) are polynomials in x, P n G .M[x]. 

Remark 1. For one- dimensional V , P n {x) are the classical orthogonal polynomials. 

Remark 2. In the construction of the present paperV can also be infinite- dimensional 
(and then M. is the Banach space of bounded linear operators on V). For simplicity 
we restrict here to matrices. 
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Relations (j3J), (jlj) are used in [5] in a different order of multiplication which is, of 
course, not essential, since it can be changed by transposition. 

The representation ([S]), ([6]) is ne and is instrumental in the deduction of the 
properties usually associated with orthogonal polynomials, §2.41 

Notation. For simplicity we denote 4- = d x . 

Proposition 1. Denote by Ak the following linear operators on M.(x) 

(5) A k = kQ'(x) + xL x + L 2 + Q(x)d x 
Then the matrices P„ defined by OD satisfy 

(6) P n = A 1 A 2 ...A n I 
forn > 1. (I denotes the identity matrix.) 

Proof. 

Calculating successively the n derivatives in (TJJ we obtain, using (j3J) and (j3J), 

d 



then 



^-[Q(x) n W(x)} = Q(x) n - 1 W(x)A n I 



^- 2 [Q(x) n W{x)} = Q(x) n - 2 W(x)A n ^A n I 



and so on. □ 

2.2. The weights W(x) defined by (J3). 

(0) The case deg Q = corresponds in the one- dimensional case to the Hermite 
polynomials, for which the operators (JSJ) are Ak = — 2x + d x . For the general matrix 
setting, besides the general properties proved for all Q, the special features of this 
case will not be considered in the present paper. 

(1) The case degQ = 1, say Q(x) = x, corresponds in the one- dimensional case 
to the Laguerre polynomials and have the representation (jSJ), © with Ak = 
—x + (k + a) + xd x . 

In the general matrix setting equation ([3]) is in this case 

(7) W'{x) = W(x) \-L 2 + U 

\ x 

If the eigenvalues of L 2 are nonresonant, in the sense that no two eigenvalues differ 
by an integer, then the general solution of ($7$ has the form W(x) = x i2 $(x) with 
analytic [9]. The behavior at infinity is exponential if no eigenvalue of L\ is 
zero. (See gH §42 for details.) 



however, a special instance was used by the author in [7]. 
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(ii) The case deg Q = 2 corresponds in dimension one to the Jacobi polynomials. 
Choosing Q(x) = x 2 — 1, the Jacobi polynomials Pjf'^ have the representation ([5]), 
(ED with 

(8) A k = (2k + a + (3)x + (a - (3) + (x 2 - l)d x 

In the general matrix setting equation ([3]) is 
(9) 

W'(x) = W(x) ( -r^—A + -^—b) with A = I {-Lt - L 2 ) , B = \(L l -L 2 ) 
\ 1 — x x + 1 / 2 2 

with general solution W(x) = (1 — x) A & + (x) = (x + l) B <&_(x) where $± are analytic 
at x = ±1 (see §4. II for details). 

(hi) An exotic type of polynomials with Q = x 2 — 1 appears naturally in a problem 
of classification of ordinary differential equations in singular regions and was studied 
i" 0- 

(iv) The commutative cases. In dimension two, if L\ and L 2 commute, it is easy 
to see that one of the following cases occurs in a suitable basis of V: 

(a) L\ and L 2 are both diagonal; 

(b) one of Li is a multiple of the identity (this case can be regarded as a particular 
instance of (a) or (c)); 

(c) one of Li is a Jordan block, and the other one, an "almost" Jordan block: 

l * = (oa)' L > = (o I) 

Note that in this last case, as a consequence of PropositionfH all the polynomials 
generated by (jl]) belong to the commutative algebra of matrices 

(10) C = {(« £); a,P^c] 

In higher dimensions the commutative cases are built up of blocks as described 
above. 

2.3. Self-adjoint weights W(x). Since orthogonality is relative to an interval J, 
then self-adjointness W(x) = W(x)* should be required on this interval (only); in 
what follows it will be assumed that x G J = [0, +oo) in the case Q = x, respectively 
xe J = [-1,1] forQ = x 2 -l. 

Proposition 2. Consider the matrix W(x) satisfying the Pearson formula (T3|). 

Assume the following nonresonance condition: (i) in the case Q(x) = x the eigen- 
values of L 2 do not differ by an integer, respectively, (ii) in the case Q(x) = x 2 — 1 
the eigenvalues of —L% — L 2 , or those of L\ — L 2 , or do not differ by an integer. 
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IfW(x) = W(x)* for x & J then W(x) reduces to scalar weights. 
The proof of Proposition^ in found in §4.31 



Proposition 3. In the two-dimensional case, dim V = 2, the conclusions of Propo- 
sition\Ehold under the more general assumption that L 2 is diagonalizable in case (i), 
respectively —L\ — L 2 , or L\ — L 2 , is diagonalizable in case (ii). 

The proof of Proposition^ in found in §4.41 

Proposition 4. Assume dim V = 2. 

The following are the only solutions of self-adjoint for x G J, not reducible to 
scalar weights. 

Define the matrices D, L, T by 

(11) 

L=(q f)> D= (o a)' T= (c d)' ^tha,(3,\,c,deR,c^0 

Let S be an invertible matrix and let L\ = S^^^LS and L 2 = S~ 1 DS . 
The following are self-adjoint solutions of (TJj).' 
(%) for Q = x: 

(12) W(x) = S*TSe xLl x L2 = e ax x x S* ( ° ° /Q , , , , , J S 

y ' y ' yc c{px + mx) + d J 

(ii) for Q = x 2 — 1 : 
(13) 

W(x) = S*TS (1 - x) A (l + x) B = (1 - x)~ x ~ a (l + x) x ~ a S* ( ° C ( . j S 



c cf(x) + d 

where f(x) = (-/? - 1) ln(l - x) + {/3 - 1) ln(l + x) 
The proof of Proposition!!] is found in §4.51 

Since all the weights f|T2|) . ([TBI have eigenvalues of opposite sign we have the 
following 

Corollary 5. //dim V = 2 any self-adjoint and positive semidefinite solution of 
with deg Q > reduces to scalar weights. 

2.4. General Properties. 

2.4.1. Commutation relations. The commutation relations (|T4^) - ([T6|) are key to the 
beautiful properties of the classical orthogonal polynomials: 

Proposition 6. For r G M.{x) the operators Ak satisfy the identities 
(14) A k (xr) = xA k r + Qr 
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(15) QA k r = Ak-i(Qr) 

(16) d x A k r = A k+ id x r + {2ak + L 1 )r 
These identities follow by an immediate calculation. 

2.4.2. Recurrence relation. 

Proposition 7. The matrix-valued polynomials P n (x) satisfy the following two-step 
recurrence relation 

(17) x P n (x) = P n+1 (x)a n + P n (x) (3 n + P n __i(x)7 n 
whith a n ,/5 n ,7„ G M. uniquely determined; in particular 

(18) a n = [Li + (2n + l)t7] -1 [L x + (2n + 2)<7p 1 [L x + (w + l)a] 
Proof of Proposition^ 

Let g G .M. The identities ( 1T4"1) . (|15l) give by iteration 

AA . . . 9) = A A • • • A1-1 At?) + • • • Ai-i (Q ?) 

(19) = ... + nAA...X_i(gg) 
Using (EJ) and ( fT9l relation ( fTTj) follows if we have 

which expanded yields an identity of quadratic polynomials in x, and by identifying 
the coefficients we obtain the following equations for a n ,/3 n ,7 n : 

(20) L x + (n + l)o- = [Li + (2n + l)<r] [Li + (2n + 2)<r] ct„ 

(21) L 2 + r = (2na + L 1 )P n 

+ [2(2n + l)<rL 2 + 2(n + 1)tLi + 2(n + l)(2n + l)ar + L X L 2 + L 2 Li] a n 

(22) - (n - 1)5 = 7 n 

+ {(nr + L 2 ) [(n + l)r + L 2 ] + 5 [2(n + l)a + Li]} a n + (nr + L 2 )/5 n 

Due to condition the system (l2"0"j) - (l2"Tj) has a unique solution a n ,/3 n ,7 n . □ 
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2.4.3. A second order equation for which P n are eigenfunctions in the commutative 
case. 

Proposition 8. If L\ and L 2 commute then P n are eigenfunctions for the operator 
A\d x and 

(23) A x d x P n =n[(n+ I) a + L x ] P n 

Proof. Using the representation ([6]) and ({16]) iteratively we obtain 

Aid x P n = Aid x AiA 2 ■ ■ ■ AJ = (2a + Li)P n + AiA 2 d x A 3 ■ ■ ■ AJ 

= ... = \(2a + Lx) + . . . + (2na + L x )\ P n 
which is (El. □ 



2.4.4. The dominant coefficient. It is clear that P n (x) are polynomials of degree at 
most n, as an immediate consequence of the representation (jSJ). The assumption (j2J) 
ensures that the degree of P n (x) is precisely n: 

Proposition 9. The coefficient of x n in P n (x) (n > 1) is 

(24) C n = (Li + 2na) [L x + (2n - l)a] . . . [L x + (n + l)a] 

which is invertible under the assumption (TJjj. 

Proof. The dominant term of P n (x) is found by retaining only the dominant terms 
in (jSJ), and using (EJ): 

P n (x) = (2ax + xL x + ax 2 d x ) . . . [(2n - 2)ax + xL x + ax 2 d x ] (2nax + xL x )+0 (x 71 ' 1 ) 

and using the fact that (XxL x + ax 2 d x )x k = (XL X + ak)x k+l the coefficient (j24"|) is 
found after a short iterative calculation. □ 

2.4.5. Completeness. The sequence {P n } n=0 12 is complete in the following sense: 

Proposition 10. For any p G M.[x] a polynomial degree n, there exist q , . . . ,q n e 
Ai so that 

n 

(25) p(x) = ^ Pk(x)q k 
and the representation [2o*\) is unique. 



Proof. The decomposition (1251) follows easily by induction on n, relying on the 
fact that the dominant coefficients (l2"4l are invertible. □ 
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2.4.6. Orthogonality. In this section the interval J is chosen so that Q(x)W(x) van- 
ishes at its endpoints. 

Proposition 11. The following quasi- orthogonality relations hold: 

(26) J P j {x)*W{x)P k (x)dx = forj<k 
and 

(27) J P j (xYW(x)*P k (x)dx = for j > k 

provided that the integrals exist. 

Relations (1261) and (1271) are an immediate consequence of the Rodrigues formula 
fll]) using integration by parts (which holds for matrix multiplication). 

2.4.7. Ladder equations in the commutative case. The ladder relations for orthogonal 
polynomials have the form 

L n P n (x) = a n P n -i(x) 

where L n are first order linear differential operators and a n are constants. These type 
of relations for orthogonal families of polynomials are found in [8] for quite general 
weight functions. In this section it is shown that in the commutative case, if the 
sequence P n satisfies a scalar type Rodrigues' formula, then the representation (j5J), 
([BJ yields easily the ladder relations. 
Indeed, consider operators of the form 

(28) L = Ax + B + CQ(x)d x , where A, B, C G C 

Note that the operators (1281) commute with A k modulo a multiple of Q(x): 
LA k = A k L + 9 k Q(x) , where 9 k = 2akC + CL X - A 

We then have 

L n Pn — L n A\A2 ■ ■ ■ A n = A\L n A2 ■ ■ ■ A n + 61QA2 ■ ■ ■ A n 

and using (}T5|) we obtain further 

= A\L n A2 ■ ■ ■ An + OiAl ■ ■ ■ An-lQ 

and iterating 

= AlA2L„A 3 ■ ..An + (61 + 62) A X ... An-lQ 

until 

= A 2 ■ -.An-l [A n L n + (&1 + . . . + 6 n ) Q ] 

The numbers A,B,C are now determined so that 

(29) A n L n + (61 + . . . + 6 n ) Q = 1 
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Expanding (|29|) we obtain a polynomial second degree in x, whose coefficients must 
vanish: 



(r + L 2 ) A + (2 na + L x ) B + n (na + a + L x ) r C = 
-7 (n — l)A + (nr + L 2 ) B + n (na + a + L x ) jC = 1 
This system has a unique solution for n > 1 under the assumption (f2l) : 

5= [a L 2 -Tna -tLx] G' 1 , C = - [2 na + L x ] G' 1 , A = -naC 



2.5. Commutative families of quasi-orthogonal polynomials. Let L\ and L 2 

belong, in a suitable basis of V, to the commutative algebra ffTOj) : 



(not both diagonal). Then P n (x) G C for all n. This is a genuine matrix setting, 
since there is no basis of V in which these polynomials are diagonal - in other words, 
the algebra C is not semisimple. 

The polynomials P n satisfy the general properties of orthogonal polynomials, as 
shown in §2.41 What is missing is orthogonality in the proper sense! 

There exists, however, a signed weight (|12|) . or ([TBI with respect to which the 
polynomials are quasi-orthogonal. 

An approach serving the function of orthogonality is the following. We will refer 
here to the case Q = x, the case Q = x 2 — 1 being similar. 

Noting that P n do not depend on the choice of the solution W(x) of (j3J) let us 
choose a solution belonging to C: taking T = / in (|T2|) (note that S = I) we obtain 



(na + a + Li) (A + naC) = 



n 



where 



G = n (4 a 7 — r 2 ) (na + Li) + 7 L 2 + a L 2 — r L x L 2 



(30) 




(31) 




Assume a.\ < and a 2 > 0. As in Proposition!!]] we have 



(32) 




Denote 



(33) 




which is obviously invertible. 
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Integration by parts gives 

roc 

(34) / P k (x) W{x) P k (x) dx = (-l) k C k M k 
Jo 

where C k is given by (|24"|) . and therefore (1341) is invertible. 

Then the coefficients q k in the decomposition (|25|) of any p £ A4[x] as an expansion 
in {P k } k >o can be found analytically as 

POO 

(35) q k = (— l) fc C k x M k x \ p(x)W(x)P k (x)dx 

Jo 

providing another justification to the use of the term "quasi-orthogonal" family. 

Note that the matrices P n {x) have classical orthogonal polynomials on the diago- 
nal. 

Generalization to higher dimensions can be made using the structure of commu- 
tative algebras which are not semisimple. 

Note that the algebra C can be re- written as the set of matrices of the form aI+/3N 
where N is a nilpotent, N 2 = 0, and a, [3 £ C. In higher dimensions one can consider 
commutative algebras consisting of elements Ylt^o a kN k where N is nilpotent with 
N n = 0, and a k 6 C. ^ ' 
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4. Appendix 

4.1. Solutions of linear systems near regular singular points - Probenius 
series. A linear differential equation in C d with a regular singularity placed, say, 
at x = has the form u'(x) = -M(x)u(x) where the matrix M(x) is holomorphic 
at x = 0. If the eigenvalues of M(0) are nonresonant, in the sense that no two 
eigenvalues differ by an integer, then the fundamental solution has the form U (x) = 
&(x)x M( - ^ with $(x) holomorphic at x = (see [9], §2.4). The matrix $ satisfies the 
differential equation $' = ^<3> (M — M(0)) which has formal power series solution 
$(x) = I + 0(x), and this series converges. 

In resonant cases the fundamental solution is also a convergent series, only it may 
also contain logarithms besides besides the (noninteger) powers (see also [TO]). 

These convergent representations of solutions at a regular singularity are called 
Frobenius series. 
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In the case (i) of §2.21 the matrix M(x) has the form M(x) = xL\ + L 2 , and, since 
the only singularities of the linear differential equation are x = and x = oo the 
function $(x) is entire. 

In the case (ii) we use M(x) = A + jrf -B to write the Frobenius series at x — 1, 
respectively M(x) = B + fr^^4 for the Frobenius series at x = — 1. The point at 
infinity is also a regular singularity: substituting x = l/£ in the equation (x 2 — l)u' = 
xLi + L 2 we obtain ^| = | |rry(£i + £,L 2 )u and the general solution has the form 
u = 1 with \1/ analytic at £ = if the eigenvalues of L\ are nonresonant. 

4.2. Local behavior of linear systems near irregular singular points. For 

systems u'(x) = (Li + ^L 2 )u the point x = oo is an irregular singularity. Assuming 
the nonresonance condition that the eigenvalues of L\ are distinct, and working in the 
basis where L\ is diagonal, there exists a transformation of the type u = (J + -S)\ 
after which the equation becomes v' = (Li + + ^g{x))w with D diagonal and g 
analytic at infinity. The fundamental system of solutions has the representation v = 
^/(l/x)x D exp(Lxx) where ^ is a formal power series in - which is Borel summable 
in generic cases, on appropriate sectors in the complex plane (see pU] for formal 
solutions, [TTJ for Borel summability in the linear case, and for nonlinear equations 
see [12] and the references therein). 

4.3. Proof of Proposition^ (i) The case Q = x. By the nonresonance assump- 
tion, in particular the eigenvalues of L 2 are distinct, therefore after a change of basis 
of V the matrix L 2 can be assumed diagonal: 

L 2 = D = diag {Ai, . . . , \ d }, (A, ^ Xj for i ^ j) 

and the equation (jHJ) defining W is 

(36) W = W (^D + Lt^j 

Using the structure of solutions at a regular singular point, see §4.11 there is a solution 
of the form Wq = x D Q(x) with $ analytic at x — 0. From (136]) we obtain that 

(37) $' + -(£>$-$£>) = $Li 

x 

which has an analytic solution $(x) = I + x$i + 0(x 2 ) with $x satisfying 

(38) $i + £>$i - $iD = Li 

Clearly any solution of (|36|) can be written as W(a;) = XWo(x) for a suitable 
constant matrix K. 

Assume that W(x) is a solution which is self-adjoint for x > 0. Then each term 
in its power expansion W(x) = Kx D (l + x$i + 0(x 2 )) must be self-adjoint. 
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For the first term we have Kx = x K* which implies KuX 4 = Kux 1 for all 
x > hence Ka, A, G JR.. Also we must have Kijx Xj = KjiX Xi for all x > and all 
z 7^ j which implies = for i ^ j, therefore K is diagonal and real. 

The auto-adjointness of the next term, Kx D $i implies that $i is diagonal and 
real by similar arguments. Finally, relation (|38|) implies that L\ is diagonal and real, 
which completes the proof. 

(ii) The case Q = x 2 — 1 follows in a similar way, by writing the series of W(x) at 
x = 1 if Li + L 2 is nonresonant, respectively at x = — 1 if L\ — L 2 is nonresonant. 

4.4. Proof of Proposition^ 

Note that if the matrices assumed diagonalizable have equal eigenvalues, then the 
proof is exactly as the proof of Proposition^ 

In the other resonant solution of (|36|) has the form Wq = x D $>(x), where now 

is a convergent series containing both x and In a: [10]. The type of logarithmic 
terms is determined by the resonant terms. 

Consider for example the case 

L 2 = D = diag{Ai, A2} with A2 — Ai = 1 
In this case W = x D $(x) is a solution if 

(39) = I + ^x n <$> n + ^2c n x n \nx E 21 

n>l n>l 

where E 2 \ is the matrix with entry 1 in the position (2, 1) and zero everywhere else. 
Substituting fl39l into (J3"7j) we obtain that $1 must satisfy 

(40) $1 + £>$i - $iD + c x E 2X = Li 

which has a unique solution when c\ equals the entry (2,1) of L\. From this point 
on the proof follows as in §4.31 

The general case Aj — Xj E Z\0 is similar, with the position of In a; in ( 1391 modified 
according to the resonance - see [TU] for a detailed description of Frobenius series. 

4.5. Proof of Proposition^. 

(i) Noting that L\ and L 2 commute then clearly f[T2"j) are solutions of ([3]). It only 
remains to show that these are the only self-adjoint solutions not reducing to the 
scalar case. 

Using Propositions^ and[3] the weights not reducing to the scalar case can only be 
found when L 2 is not diagonalizable, therefore when L 2 = S~ l DS with D as in pip , 
for some S. 

The line of the proof is similar to the one of §4.31 (only with a different outcome!). 
Details are provided below, for completeness. 
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Denote L = SLtS' 1 . Equation © becomes W' = W (±D + L) where W = WS, 
therefore W = Ke D $(x) with $ satisfying (J37J), = I + x$i + 0(x 2 ). We have 
W = Ke D $S. 

Denote T = S~ l K (a constant matrix); the relation W = W* is equivalent to the 
self-adjointness of Te D $ for x > 0, which implies self-adjoint ness of all the terms in 
the series, in particular of the first term Tx D , and of the second term Tx D $i. 

We have 

d a ( 1 In x 

X = X ' 



1 

and a simple calculations shows that Tx D is self-adjoint only for AeM and T as in 

(HU). 

Then Tx D <&\ is self-adjoint only for matrices of the form 

* = {o I )• a - i3eR 

Since <3>i must solve (|38|) (with Li = L), it follows that L = $i, which completes 
the proof of Proposition|4] in the case Q = x. 

(ii) The case Q = x 2 — 1 is similar. □ 
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